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Abstract. We briefly outline the physics underlying the mechanical and thermal 
f^ • properties of brown dwarfs, which characterizes their interiors and their atmospheres. 

5^ , We mention the most recent improvements realyzed in the theory of brown dwarfs and 

the connection with experimental and observational tests of this theory. 
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•^ ! 1. Introduction 

o 

Cs| ■ A general outline of the basic physics entering the structure and the evolution of brown 

Qv ' dwarfs (BD) can be found in the previous reviews of Stevenson (1991) and Burrows 

0^ ■ & Liebert (1993). Important innovations have occured in the field since then, one 

of the least negligible being the discovery of bona- fide brown dwarfs (Rebolo et al., 
1995; Oppenheimer et al., 1995). An increasing number of these objects have now 
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i-i • been discovered either as companions of stars, as members of young clusters or as free 

^ ■ floating objects in the Galactic field (Ruiz et al., 1997; Delfosse et al., 1997). On the 

other hand the theory has improved substantially within the past few years and can now 
L^ ■ be confronted directly to observations and even to laboratory experiments, as will be 

Vh ' shown below. It is thus important to reconsider the previous reviews in the light of these 

observational and theoretical progress and to update our knowledge of the structure and 
the evolution of BDs. This is the aim of the present review. 

BDs are objects not massive enough to sustain hydrogen burning in their core 
and thus to reach thermal equilibrium, defined as L = Lnuc where Lnuc = Jo edm 
is the nuclear luminosity and e is the nuclear reaction rate per unit mass. This 
hydrogen burning minimum mass Mhbmm depends on the internal composition of 
the object, in particular the abundance (by mass) of helium (Y) and heavier elements 
(Z). For abundances characteristic of the solar composition, typical of the Galactic 
disk population, Yq=Q.27 and ^©=0.02, this minimum mass is Mhbmm ~ 0.072 M©, 
whereas for compositions characteristic of the Galactic halo (Y=0.25; Z ~ 10~^ x Zq) 
Mhbmm ~ 0.083 M0 (Chabrier & Baraffe, 1997; Baraffe et al., 1997) and for the zero- 
metalhcity limit {Z = 0), Mhbmm ~ 0.09 M© (Saumon et al ., 1994). 

The minimum mass for BDs is presently undetermined and masses as small as a 
Jupiter mass {IO^^Mq) are not excluded in principle. The divided line between BDs and 
giant planets (GPs) is still unclear and stems essentially from their formation processes : 
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hydrodynamic collapse of an interstellar molecular cloud for BDs, like for stars, accretion 
of heavy elements in a protostellar disk for the formation of planetesimals which 
eventually become dense enough to capture hydrogen and helium and form gazeous 
planets. The accretion scenario rather than the collapse scenario for the formation of 
planets is supported by the distinctly supersolar average abundance of heavy elements 
in Jupiter and Saturn, although there is only indirect evidence for the presence of the 
central rocky core through the modeling of the centrifugal moments. The border line 
between these two scenarios is presently unknown and involves most likely complex 
dynamical and non-linear effects. Extra-solar planets with masses as large as ~ 40 Mj 
and BDs with masses as low as ~ 30 Mj have now been discovered. Except for this 
formation process and for the presence of a central rock/ice core, the physics and the 
observational signature of BDs and GPs is very similar. Since a complete review is 
devoted to GPs (Stevenson, this issue), the present one will be devoted to BDs. In the 
present survey, I will focus on the most recent improvements realyzed in the physics 
of the interior and the atmosphere of BDs. I will also mention the physics underlying 
the so-called Lithium-test, which provides a powerful independent determination of the 
substellar nature and the age of a putative BD. The aim of the present review is not to 
present detailed calculations (which can be found in the various mentioned references) 
but rather to catch the underlying physics entering the structure and the evolution of 
BDs. 

2. Interior of brown dwarfs. The hydrogen equation of state 

Central conditions for massive BDs are typically Tc< 10^ K and p^ ~ 10^-10^ 
g.cm^^. Under these conditions, the average ion electrostatic energy {Ze'^)/a, where 
a = {^-^Y^^ is the mean interionic distance, is several times the average kinetic 
energy kT, characterizing a strongly coupled ionic plasma with a coupling parameter 
Ti = {ZeY/akT > 1. The temperature is of the order of the electron Fermi temperature 
kTp and the average inter-electronic distance ag is of the order of both the Bohr radius 
tte ~ ao and the Thomas- Fermi screening length Og ~ otf- We thus have to deal with a 
partially degenerate, strongly correlated, polarizable electron fluid. The temperature in 
the envelope is kT < 1 Ryd, so we expect electronic and atomic recombination to take 
place. At last the electron average binding energy is or the order of the Fermi energy 
Ze'^/aQ ~ ep so that jjress'ure-ionization is taking place along the internal density profile. 
Recently laser-driven shock-wave experiments have been conducted at Livermore 
(Da Silva et al., 1997; Collins et al , 1998) which probe directly the thermodynamic 
properties of dense hydrogen under conditions characteristic of BDs and GPs. The 
relevance of the Livermore experiments for the interior of these objects can be grasped 
from Figure 1. About a decade ago, Saumon & Chabrier (1991, 1992; SC) have 
developed a free energy model aimed at describing the thermodynamic properties of 
strongly interacting H2 molecules, H atoms, H+ protons and electrons under such 
astrophysical conditions. The abundances of each species derive from the free energy 
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Figure 1. Phase diagram of hydrogen in T — p. The coexistence curve of the 
plasma phase transition (PPT) appears at the left of center as a black solid line which 
ends at the critical point. Curves of constant ionic (P) and electronic (vs) plasma 
coupling parameter and electron degeneracy parameter 9 = T/Tp are shown. Regions 
dominated by molecules, atoms and ionized H are labeled and delimited by a curve 
which corresponds to 50% dissociation and ionization. The calculated shock Hugoniot 
corresponding to the experiments is shown by the long-dash curve across the diagram. 
Finally, the thin dotted lines show the internal structure profiles of several astrophysical 
bodies (from left to right), Jupiter, a 7 Ga brown dwarf of 0.055 M0 and a 0.1 Mq 
star. The EOS model is invalid in the hashed region. 
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This model relies on the so-called chemical picture, which assumes that the species 
remain distinct even at high-density. This requires the knowledge of the interparticle 
potentials (Ph2H2^(Ph2H and (pHH- Since N-body effects strongly modify the interaction 
between particles at high-density, "effective" pair-potentials can be derived from the 
experimental hugoniots. These effective pair-potentials mimic the softening of the 
interaction due to the surrounding particles and thus retain some density-dependence in 
the characteristic interactions between the main species. At the time the SC formalism 
was derived, the only available shock-wave experiments were the ones by Nellis et 
al. (1983). Since molecular dissociation was negligible under these conditions, only 
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an effective potential between H2 molecules could be derived. The H2-H and H-H 
potentials were taken as ab-initio potentials. The more recent experiments by Weir 
et al. (1996) reach higher pressures and substantial molecular dissociation is inferred 
from these experiments {X^ > 20%). This allows us to derive effective potentials for 
(PH2H2 ^^^ 4>H2H as well and thus to update the SC model. 

The recent laser-driven experiments have shown that the agreement between the 
predictions of the SC model and the data is excellent. In particular the strong 
compression factor arising from hydrogen pressure-dissociation and ionization observed 
in the experiment {p/pi ~ 5.8) agrees well with the predicted theoretical value (Saumon 
et al., 1998). The compression is slightly underestimated in the theory and starts at 
a slightly too large pressure. This reflects the underestimated degree of dissociation 
in the model, which stems from the too repulsive (ab-initio) 0^2^ potential at the 
time the SC model was elaborated. This shortcomings is resolved when including the 
afore-mentioned new effective H2-H2 and H2-H potentials. Eventually full ionization 
is reached at very high pressure (P ~ 10 Mbar), characterized by the asymptotic 
compression factor Pf/pi = 4 for a monoatomic fully dissociated proton fluid. These 
results show that, although this "chemical" model certainly does not pretend to give an 
exact, complete description of all the interactions in the high-pressure strongly correlated 
fluid, it very likely retains the main physics underlying the phenomenon of pressure- 
dissociation/ionization. 

It is worthnoting that a similar strong compression factor is obtained also with the 
so-called fugacity expansion scheme, in principle exact in the strongly dissociated regime 
(Rogers and Young 1997), although this scheme fails at lower density when substantial 
recombination occurs. 

One of the most striking features of the SC theory is the prediction of a first- 
order so-called plasma phase transition (PPT) between a molecular state and a plasma 
state for the pressure-ionization of hydrogen, similar to the one predicted originally by 
Wigner & Huntington (1935). However it is important to stress that the PPT in the 
SC model arises from first-principle thermodynamical instability of the one single free 
energy model {{dP/dp)T < 0) and not from the comparison between two different free 
energy models. The new SC EOS, incorporating the new potentials, still predicts a 
PPT, although with a critical point slightly cooler than predicted previously, namely 
Tc = 14600 K, Pc = 0.73 Mbar (Saumon et al., 1998). In order to really nail down the 
existence of the PPT, we have calculated a second-shock Hugoniot reflected from the 
principal one, which should be realizable in a near future (Saumon et al., 1998). Such 
an experiment should conflrm or rule out deflnitely the presence of the PPT. 

The main question about the PPT is : if it exists, what is its nature ? This question 
has been addressed to some extend in Saumon & Chabrier (1992). If the PPT exists, 
it stems very likely from the large difference between a molecular state characterized 
by a strongly repulsive potential and a plasma state characterized by a soft Yukawa- 
like potential. Given the large difference between these two potentials, and thus the 
respective available phase spaces, we can expect a discontinuity in the interaction energy 
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and thus an abrupt change in the two-particle distribution function. This behaviour is 
observed in recent path-integral Monte-Carlo simulations (Magro et al., 1996; Ceperley, 
this issue). In terms of ground state energies, this translates into the large energy 
barrier between the ground state energy of an H2-like system (H2 or Hg") and an H"'"-like 
system. In terms of correlation lengths that characterize the many-body effects, the 
system will collapse from a dense molecular phase characterized by a length Xh2 ~ a 
few Qq into a plasma phase characterized by a length \h+ « \h2- The underlying 
critical quantity will be the electron correlation length, with a critical percolation 
from a "bound-electron" -like value to a "free-electron" -like value. In this sense the 
PPT resembles the metal-insulator transition in metals associated with the liquid- 
vapor transition (Hensel, this issue), leading eventually to a polarization catastrophe 
(Goldstein & Ashcroft, 1985). The effect is likely to be more dramatic for hydrogen 
because of the absence of core electrons. 

In this sense, the conductivity measurements of dense fluid hydrogen by Weir 
et al. (1996) do not rule out the PPT. The conductivity exhibits a plateau with 
a ~ 2000 {Q.cm)^^ up to the highest pressure reached, P ~ 1.8 Mbar. This is still 
orders of magnitude smaller than the conductivity characteristic of a fully dissociated 
plasma phase, a ~ 10^ {Q.cm)~^ (Stevenson & Ashcroft, 1974) and is consistent with 
conduction being due to delocalized electrons from H2. This does not preclude a 
structural transition like the PPT at higher pressures. 

If the PPT exists it can have important consequences for BDs and GPs. The 
interior of these objects is essentially isentropic. Since the signature of a first-order 
transition is a density and entropy discontinuity, integration along the internal adiabat 
from the observed outer conditions yields different central conditions with and without 
PPT (Chabrier et al., 1992). In principle the signature of the PPT in the interior of 
GPs like Jupiter and Saturn could be observed from the analysis of p-niode oscillations 
(Marley, 1994; Gudkova et al., 1995). However this requires very accurate observations 
of high-degree modes, a difficult observational task. The PPT bears also important 
consequences on the evolution of these objects. Since by definition BDs and GPs 
do not sustain hydrogen burning, application of the first and second principles of 
thermodynamics yields the following equation for their evolution : 

dt Jo ^ p^ dt' Jo dt ^ ^ 

where L is the luminosity, u and s the specific internal energy and entropy, 
respectively. If the PPT exists, an additional term, namely the latent heat of the 
phase transition, must be added to the previous equation : 

L' = L+ f T^dm (3) 

J Am dt 

This effect was first pointed out by Stevenson & Salpeter (1977) and examined in 
detail by Saumon et al. (1992). 
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3. The atmosphere of brown dwarfs 

3.1. Spectral distribution 

The photosphere is defined as the location where the photon mean free path is of the 
order of the mean interparticle distance, i.e. It, ~ l/{Rp) ~ a oc p~^^^, where R ~ Icrn^ / g 
is the mean absorption coefficient (opacity). This equahty yields /j, ~ a ~ 1 cm. In terms 
of the dimensionless optical depth r = z/lu, where z is the depth of the atmosphere, 
equilibrium between internal and gravitational pressure yieds : 

_dP 

dr = — (p/t) dz = K, — , (4) 

where g = GM/R^ is the surface gravity. For BDs, M <Q.IMq, R ~ 0.1 i?©, 
(?<10 X gQ. This yields Pph ~ g/R ~ 10 bar at the photosphere, and pph ~ 
10^^ — 10~"^g.cm"^. Collision effects are signiffcant under these conditions. Therefore 
thermodynamic equilibrium can be safely assumed near the photosphere. The bad news 
is that collision effects can induce dipoles between molecules, e.g. H2 or He-H2, which 
otherwise would have only quadrupolar transitions. This so-called collision-induced 
absorption (CIA) between roto-vibrational states {y — > v') of e.g. 2 H2 molecules (1 and 
2) can be written in terms of the 2-body absorption (see e.g. Borysow et al., 1985): 

I^H2H2 — ^Vi,v[^V2,v'2'^H2H2 l^'-* ) 

27r2 
= <^«^(1 - e-""/'^)S.„.iS.„„.^-'<'-''^Hc.,T) (5) 

where uj = 2tiv is the angular frequency, nH2 is the number density of hydrogen 
molecules and g'"^''"^''"'^''"^2{lJJ,T) is the spectral function. We note the dependence 
on the square of the number abundance. As temperature decreases below ~ 4000 
K, an increasing number of hydrogen molecules form and thus H2 CIA-absorption 
becomes overwhelmingly important, a feature shared with giant planet and white dwarf 
atmospheres. Since the CIA absorption of II2 under the conditions of interest for BDs 
and GPs takes place around 2.2 yum, energy conservation leads to a redistribution of the 
emergent radiative ffux toward shorter wavelengths (Saumon et al., 1994; Baraffe et al., 
1997). 

The effective temperature is defined as the integral of the Eddington flux over the 
frequency spectrum : 

T^f, = a-' J H^du, (6) 

where a = 5.67 x 10^^ erg.cm^.K'^.s^^is the Stefan-Boltzman constant. BDs are 
characterized by effective temperatures Teg < 2000 K. At these temperatures, numerous 
molecules like e.g. H2, H2O, TiO, VO, ... are stable and are the major absorbers of 
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photons. These strongly frequency-dependent opacity sources yield a strong departure 
from a black-body energy distribution (see e.g. Figure 5 of Allard et al. (1997)). An 
updated detailed review of the physics of the atmosphere of low-mass stars and BDs 
can be found in Allard et al. (1997). 

Below T ~ 1800 K, carbon monoxide CO is predicted to dissociate and to form 
methane, CH4, as observed in Jupiter. This prediction has been confirmed by the 
discovery and the spectroscopic observation of Gliese229B. The presence of methane in 
its spectrum assessed unambigously its sub-stellar nature. Consistent synthetic spectra 
and evolutionary calculations done both by the Lyon group (Allard et al., 1996) and 
the Tucson group (Marley et al., 1996) yielded the mass determination of the object 
between ~20 and 50 Mj, the undermination in the mass reflecting the undetermination 
in the age of the system. 

At last below ~ 2000 K complex compounds (grains, also called "clouds" by 
planetologists) condensate in the atmosphere (see e.g. Lunine et al., 1986; Tsuji et al., 
1996). These grains will affect the atmosphere in different ways. They first modify the 
EOS itself and thus the atmospheric temperature/density-profile, and they also strongly 
affect the atmospheric opacity and thus the emergent radiation spectrum. At last they 
will produce an increase of the temperature in the uppermost layers of the atmosphere, 
the so-called backwarming (or greenhouse) effect, destroying otherwise stable polyatomic 
species. The condensation of the grains in a BD atmosphere is illustrated in Figure 
2. Spectroscopic observations of different BDs at various effective temperatures show 
evidence for an even more complicated problem, namely grain diffusion (settling) in the 
atmosphere. 

3.2. Energy transport 

The radiative transport equation reads : 

A d , VT 
Frad = ^r——[(yT ) oc ^ (7) 

6Kp dr K 

for the radiative fiux, while convective transport can be writen as : 



Fconv OC {pVconv) X {Cp5T) (8) 

where Vconv is the convection velocity, typically a fraction of the speed of sound, Cp is 
the matter specific heat at constant pressure and 5T is the energy difference between the 
convective eddy and the surrounding ambient medium. As the temperature decreases 
below ~ 5000 K, which corresponds to a mass m < 0.6 Mq, H atoms recombine, nH2 
increases, and so does R through H2 CIA-absorption (see above). The opacity increases 
by several orders of magnitude over a factor 2 in temperature. On the other hand, the 
presence of molecules increases the number of internal degrees of freedom (vibration. 
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Figure 2. Relative abundances of gas phase (full lines) and crystallized species 
(dotted lines) across a Toff = 1800 K brown dwarf model atmosphere (after F. AUard). 



rotation, electronic levels) and thus Cp. These combined effects strongly favor the onset 
of convection in the optically-thin (r < 1) atmospheric layers. This can be shown easily 
from a stability (Schwarzchild) criterion analysis. Flux conservation thus reads : 



V{Frad + F^ 







(9) 



i.e. no radiative equilibrium. The evolution of low-mass objects (low-mass stars, 
BDs, GPs) thus requires to solve the complete set of the transfer equations and to 
use consistent boundary conditions between the atmosphere and the interior structure 
profiles (Chabrier & Baraffe, 1997; Baraffe et al., 1995, 1997, 1998; Burrows et al., 1997). 



4. Screening factors and the Lithium-test 

Since a BD, by definition, never reaches thermal equilibrium (L ~ TdS/dt), age is an 
extra degree of freedom, yielding an undetermination in the mass and/or age of an object 
for a given observed luminosity and/or temperature. An independent age-indicator is 
thus needed. The presence of Lithium in the atmosphere of a cool object provides such 
an indication. The signature of Lithium absorption as a diagnostic for the sub-stellar 
nature of an object was first pointed out by Rebolo et al. (1992) while the measure of 
Lithium-depletion as an age tracer was first used by Basri et al. (1996). 

The physics underlying the Lithium-test roots in dense plasma physics and in the 
calculations of the so-called nuclear screening factors for the nuclear reaction rate. 
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Figure 3. Central temperature as a function of age for 3 different masses respectively 
above, at the limit and below the hydrogen-burning minimum mass. Th and T^i 
indicate the Hydrogen and Lithium burning temperatures, respectively. 



Primordial '^Li is detroyed through the nuclear reaction "^Li + p -^ 2'^He. The 
reaction rate R^ (in cni~^s~^) in the vacuum is given by the usual Gamow theory 
Rq oc e"^^"/*^^ where eo corresponds to the Gamow-peak energy for non-resonant 
reactions, which corresponds to the maximum probability for the reaction. However, 
as mentioned above, non-ideal effects dominate in the interior of BDs and lead to 
polarization effects in the plasma. These polarization effects due to the surrounding 
particles yield an enhancement of the reaction rate, as first recognized by Schatzman 
(1948) and Salpeter (1954). The distribution of particles in the plasma reads : 



n[r) 



ne 



-Ze<p(r)/kT 



(10) 



with 



Ze 



+ ■?/'(r) 



:ii) 



where %l){r) is the induced mean field potential due to the polarization of the surrounding 
particles. This induced potential lowers the Coulomb barrier between the fusing particles 
and thus yields an enhanced rate in the plasma R = E x Rq where 



E = limr-,0 {gi2ir)exp{ \^ )} 

rkl 



(12) 



is the enhancement (screening) factor and (712 (^) the pair-distribution function. 

Under BD conditions, not only ionic screening must be included but also electron 
screening, i.e. E = Ei x E^. Both effects are of the same order {Ei ~ i^g ~ a few) 
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and must be included in the calculations for a correct estimate of the Lithium-depletion 
factor [Li]o/[Li], where [Li]o = 10~^ denotes the primordial Lithium-abundance. This 
yields a Lithium-burning minimum mass rriLi ~ 0.06 M© (Chabrier & Baraffe, 1997), 
below the hydrogen-burning minimum mass, as illustrated in Figure 3. After the common 
primordial deuterium burning phase, which lasts ~ 10^ — 10^ yr, the central temperature 
evolves differently, depending on the mass of the object. Note the strong age dependence 
of the Lithium-test : young stars with an age t < 10^ yr (depending on the mass) will 
exhibit Lithium, whereas massive brown dwarfs within the mass range [0.06-0.07 Mq] 
older than ~ 10^ yr will have burned Lithium. The measure of Lithium depletion in the 
atmosphere of low-mass objects, inferred from the width of the Lil line at 6708 A, as an 
age indicator, is illustrated in Figure 4. This figure displays the evolution of a 0.075 Mq 
object, the H-burning limit for solar-abundances, in the I-band magnitude. The left and 
right diagonal solid lines correspond to 50% Li-depletion ([Li]o/[Li] = 1/2) and 99% Li- 
depletion, respectively. Thus, for say 120 Myr, the inferred age of the Pleiades cluster, 
objects brighter than Mj ~ 12.2 will lie on the right-hand side of the 99%-depletion line 
and thus are predicted to show no Lithium in their atmosphere and to be H-burning 
stars {m > 0.075 M©), whereas objects fainter than this magnitude will all show some 
Lithium and all be brown dwarfs {m < 0.075 Mq for this age), with objects fainter 
than Mj ~ 12.6 predicted to have retained more than half their primordial Lithium- 
abundance. The horizontal lines show the observed magnitudes of 4 different objects in 
the Pleiades, with available high-resolution spectra. All four confirm the theory, with no 
Lithium observed for PLIO, about 50% depletion for PL 13 and negligible or no depletion 
for Roql3 and Teidel. Different isochrones for different masses can be superimposed 
on the same diagram and analyzed similarly. This illustrates convincingly the powerful 
diagnostic of Lithium as a mass and age indicator for low-mass stars and brown dwarfs. 

5. Conclusion 

As a conclusion to this review, I will list a series of "homework problems" related to 
BDs which illustrate the major problems to be addressed in this field in a near future 
and which correspond to different domains of physics or astronomy. This list is certainly 
not exhaustive. 

• Dense matter physics: As we have seen, BD interiors can now be tested 
directly in laboratories and the EOS of these objects can be probed by high-pressure 
experiments. More experiments are needed in the complex regime of H-pressure 
dissociation/ionization with several unanswered questions. Does the PPT really exist ? 
Does it survive when 10% helium particles are present ? How does pressure-ionization 
of H affect the dynamo process in BD and GP interiors ? 

• Star formation process: Jeans stability analysis yields a minimum mass mmin ~ 
0.01 Mq, definitely in the BD domain (Silk, 1977). Is this mass the BD minimum mass ? 
Conversely what is the maximum mass for planet formation ? Does the Jeans criterion 
really apply for the formation of star-like objects ? What is the BD mass function in 
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Figure 4. Evolution of tlic absolute magnitude M/ as a function of age. The dashed 
line corresponds to the hydrogen-burning minimum mass, whereas the diagonal solid 
lines correspond to the 50% et 99% Lithium-depletion limit. The horizontal solid lines 
indicate the observed magnitudes of different low-mass objects. 



the Galaxy ? 

• Evolution: The evolution of BDs is not hampered by any adjustable parameter, 
like for example in the treatment of convection for more massive stars which develop 
an inner radiative core. The theory of BDs, and the comparison with observation, thus 
reflects the validity of the very physics entering the theory, both in the atmosphere and 
in the interior. This theory can be tested directly now by photometric and spectroscopic 
observations and must address new problems like e.g. the diffusion process of grains in 
the atmosphere or the magnetic field generation in active BDs. Conversely, the theory 
is now reliable enough to provide useful guidance for future observations. 

• Galactic implication: The mass-to-light ratio for BD, (M/L)^£) > 10^(M/L)q, 
make BDs very promizing candidates to explain at least the baryonic missing mass. 
Even though present estimates of their contribution to the Galactic disk and halo 
mass seem to exclude this possibility (Chabrier & Mera, 1997; Mera et al., 1998), the 
determination of their exact number- and mass-density in the Galaxy remains to be 
determined accurately. Ongoing microlensing experiments sensitive to hours and day 
event durations and ongoing wide field infrared projects (e.g. DENIS, 2MASS) will 
certainly help nailing down this issue. 

BDs thus present a wide variety of interest from basic physics to Galactic 
implications and should remain a very active field. 
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